Ghost spins and novel quantum critical behavior in a spin chain 

with local bond-deformation 



Jianhui Dai a,fe Yupeng Wang c,d U. Eckern c 
a Abdus Salam International Center for Theoretical Physics, Trieste 34100, Italy 
b Zhejiang Institute of Modern Physics, Zhejiang University, Hangzhou 310027, China 
c Institut fur Physik, Universitat Augsburg, D-86135, Augsburg, Germany 
d Cryogenic Laboratory, Chinese Academy of Sciences, Beijing 100080, China 



Abstract 

We study the boundary impurity-induced critical behavior in an integrable 
5C/(2)-invariant model consisting of an open Heisenberg chain of arbitrary 
spin-S 1 (Takhatajian-Babujian model) interacting with an impurity of spin S' 
located at one of the boundaries. For S = 1/2 or S' = 1/2, the impurity inter- 
action has a very simple form J Si ■ S' which describes the deformed boundary 
bond between the impurity S' and the first bulk spin Si with an arbitrary 
strength J. With a weak coupling < J < Jq/[(S + S') 2 — 1/4], the impu- 
rity is completely compensated, undercompensated, and overcompensated for 
S = S' , S > S' and S < S' as in the usual Kondo problem. While for strong 
coupling J > Jo/[(S + S') 2 — 1/4], the impurity spin is split into two ghost 
spins. Their cooperative effect leads to a variety of new critical behaviors 

with different values of |S" — S\. 
75.10.Jm, 75.20.Hr, 72.15.Qm, 05.70.Jk 
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I. INTRODUCTION 



on 



Quantum fluctuations induced by an impurity coupled to the one-dimensional Tomonaga 
Luttinger liquid (TLL) play essential role in understanding the low temperature behav- 
ior of quasi- ID systems, such as quantum wireso, fractional quantum Hall effect!, carb 
nanotubesil or quasi- ID organic conductorsS. The problem of an impurity spin S' coupled 
with both of the neighboring sites of quantum chains was studied by a class of integrable 
SU (2)-invariant modelsiil. For the Heisenberg chain with ferromagnetic coupling, the impu- 
rity is locked into the critical behavior of the lattice, i.e., at low temperatures the specific heat 
is proportional to T 1 / 2 and the susceptibility diverges as T~ 2 with logarithmic corrections^. 
For a chain with antiferromagnetic coupling, the impurity spin is compensated by bulk spins 
with three different situations similar to the multichannel Kondo problem!: for 5" = S, it is 
the complete compensation and the impurity just corresponds to one more site in the chain; 
for 5" > S, the partial compensation with Schottky anomaly when an external magnetic 
field H is applied; and for S' < S, the overcompensation which gives rise to quantum critical 
behavior!. 

Meanwhile the effects of an impurity embedded in ID TLL were extensively discussed 
recently. By renormalization group(RG) techniques, bosonization methods and boundary 
conformal field theories, many interesting results have been obtained, showing the unusual 
properties of TLL in the presence of a local potential barrier or a magnetic impurityiHil. 
Generally speaking, these new findings indicate that the quantum impurity models renor- 
malize to critical points corresponding to conformally invariant boundary conditions^Hii. 
The impurity-bulk coupling strength J flows either to infinity when the impurity is screened, 
or to finite as if it is overscreened, no matter what the sign of J is initially. In particular, 
numerical studies of the finite size spectrum support the picture that the fixed point corre- 
sponds to a chain disconnected at the impurity site for repulsive interaction^. However , the 
low temperature impurity behavior described by previous Bethe ansatz integrable models 
do not correspond to the stable critical points mentioned above. For instance, at the critical 
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point the spin-5 impurity coupled to spin-1/2 antiferromagnetic Heisenberg chain has the 
effective screened spin S e ff = S — 1/2 rather than 5 — 1, despite the fact that it couples to 
the two neighboring 1/2-spin'si. In this respect, the critical point described by the impurity 
model is unstable, owing to the fact that these models have a fixed impurity coupling, a 
"fine tuned" impurity interaction term, and no backward scattering^. It is recalled that 
backward scattering is one of the essence of quantum impurity problem in ID TLLi0. From 
the point of view of RG, electrons or spin waves moving in one -dimensional space will be 
largely scattered back by the impurity, while the tunneling effect could be perturbed00. 
In the fixed point limit, they are completely scattered back after a phase shift due to in- 
teraction with the impurity, as long as transmission being plausibly neglected at sufficiently 
low temperature. Based on the fixed point observed by RG, some open boundary problems 
with the impurity located at the boundary were consideredHB. By comparing the open 
boundary models and the periodic models, Zvyagin and co- authorUS argued that the im- 
purities in the integrable models could be generic. However, we remark either a transparent 
impurity or a boundary impurity in the integrable models considered earlier can not reveal 
a full picture of a generic impurity in the bulk since only one channel host is included in 
these models (for the transparent impurity, only forward scattering while for the boundary 
impurity, only one half-chain). To make the problem to be generic, one should consider two 
half-chains interacting symmetrically with the impurity^!. At low energy scales, the problem 
is effectively a two-channel one EE as long as the interaction in the bulk is repulsive. Gen- 
erally, tunneling through the impurity may exist. This causes hybridization, splitting and 
anisotropy of the two channels. However, the tunneling matrix is negligibly small comparing 
to the Kondo coupling[i-!l or the impurity potential! at low energy scales from the RG point 
of view. Therefore, these effects are not very harmful to the two-channel behavior@. 

In this paper, we consider the problem of an impurity spin interacting symmetrically 
with two half spin-chains. By mapping the problem to an open boundary one, we solve the 
related integrable model via Algebraic Bethe ansatz (ABA). The structure of the present 
paper is the following: In the subsequent section, we construct the model and derive the 



Bethe ansatz equation. The ground state properties and the boundary bound states will be 
discussed in sect. III. In sect. IV, we discuss the thermodynamics of the open boundary as 
well as the impurity. It is found that the open boundary behaves as an overscreened spin 
and the impurity itself, may show different quantum critical behaviors, depending on the 
coupling constant J. Sect.V is attributed to the concluding remarks. 



II. THE MODEL AND ITS BETHE ANSATZ 

Let us start with the following Hamiltoniani! 

iV-1 

H = J J2[Sn- S n+1 + £_„ • 5_ n _ x ] + JS> ■ {St + (1) 

n=l 

— * — t 

where Jo > (antiferromagnetic coupling); S n is the spin- 1/2 operator, S' is the impurity 
spin operator. As long as J ^ J or S' ^ 1/2, back scattering off the impurity dominates over 
the forward scattering and the problem is effectively a two-channel one. No matter J > 
(antiferromagnetic) or J < (ferromagnetic), the second term in (1) causes ferromagnetic 
correlation between Si and That means S\ and S"_i favors to form a spin-1 composite at 
low temperatures. In fact, there is a quasi-long-distance ferromagnetic correlation between 
S n and S_ n 

< S n ■ S. n >~ (2n)- e , (2) 

where the exponents 9 varies from 1 to 4, which can be derived from the boundary conformal 
filed theoryll!. Therefore, there is a tendency for S n and S- n to form a spin-1 composite 
and the problem could be mapped to a spin-1 chain coupled with a boundary impurity. We 
remark the composites far away from the impurity are not very tight and may lose sense as 
can be seen in (2), but this does not change very much the impurity behavior since the bulk- 
impurity interaction is local and very short-ranged. Hence, a spin-5* chain with a boundary 
impurity keeps the main feature of an impurity in a spin-S72 chain. Caution should be 
taken when we construct the spin-S* open chain to avoid the Haldane gap since the original 
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Hamiltonian (1) is gapless. The formation of the high-spin composites is very similar to 
that in a multi-channel Kondo problem! but with a different mechanism. Of course, the 
boundary impurity-induced critical behavior is quite generic. 

Based on the above discussion, we study the low temperature behavior induced by an 
magnetic impurity S' coupled to an open antiferromagnetic Heisenberg chain of spin 5 by use 
of Bethe ansatz exact solutions. It is well known that integrable generalization of isotropic 
S = 1/2 spin chain to arbitrary spin 5 leads to the Hamiltonianil 

AT-l 

Hs — Jo X] Qzs{Sj ■ Sj + i), (3) 
where Q2s(x) is a polynomial of degree 25 of SU{2) invariant quantities x = Sj ■ Sj + \ 

25 j i 25 _ 

Qis(x) = UE- k ) n (4) 

j=l fc=l ^ l^j,l=0 3 x ' 

with x n = |n(n + 1) — 5 (5 + 1), n = 0, 1, • • • , 25. One recovers Hy 2 = Jo X/^Ll 1 Sj • Sj+%, 
Hi = Jo Ysf^ilSj ■ Sj+i — (Sj ■ Sj + i) 2 ] as the usual Spin-1/2 Heisenberg model and the 5 = 1 
Takhatajian-Babujian model respectively (up to an irrelevant constant). The construction 
of the model is based on the vertex weight operators R(X), represented by matrices acting 
on the tensor product spaces V\ ® Vi of two spins Si, S 2 , with a parameter A identified as 
spin rapidity. Its explicit form is 

5* U (A) = -Enrn^> ( 5 ) 

1=0 k=0 A + K 

where Pi is the projector selecting the states with total spin I in the tensor product of the two 
spins involved, Pi(x) = Y\n%, n =o x~-x n ' ^ wm S ^° ^ ne Yang-Baxter equations satisfied by the 
_R-matrix, we have the relationship Hs oc 4- lni(A)|,\=o> with t(X) being the transfer-matrix 
defined by 

t(X) = tr A T(X) = tr A { s R AN (X) ■ ■ - s R A1 (X)}. (6) 

Here, the trace is taken in the auxiliary spin space Va (dimVA = 25 + 1) introduced to help 
us track of the proliferating spin indices. Because [t(X),t(fi)] = 0, VA, fi, Hs is integrable 
under the periodic boundary condition. 



Now, we put a magnetic impurity S' at one of the end of an open spin-S Heisenberg 
chain, by considering the following integrable Hamiltonian 

H = H s + H imp , (7) 

S+S' l u S+S> 

H imp = J E ( E t^— 2 ) II (8) 

l=\S-S'\+l k=\S-S'\+l C n^l,n=\S-S'\ ^ ^ n 

where?/ = Sx-S', yi = ~[l(l + l)—S(S+l) — S'(S'+l)]; cis an arbitrary parameter describing 
the strength of the bulk-impurity interaction. The impurity is assumed to be sited at the 
left hand end of the chain, say the site j = 0, while its neighboring spin is Si. For S = 1/2, 
the model is reduced to that considered in Ref.[21]. Interestingly when S = 1/2 or S' = 1/2, 
the interaction term takes the simple form 

Himp = J Si ■ S', (9) 

with coupling constant J = Jq/[(S' + S) 2 — c 2 ], which can range from negative infinity 
to positive one, and meet all the physical situations. So at least in these two cases, the 
Hamiltonian could be expected to properly describe the boundary bond effect in some real 
quasi-lD materials at very low temperature, such as the possible bond impurity S' = 1/2 
in S = 1 Heisenberg antiferromagnet TMNIN0. 

To show the integrability of the Hamiltonian (7), let us first notice that the impurity 
term H imp could be more conveniently treated as the boundary operator, similar to the usual 
open boundary problem with an external field applied to the end. In addition to the Yang- 
Baxter equation (YBE) as the integrable condition of the bulk, there are some new consistent 
constraints (often called the reflect YBE ) for the same model to be integrable under the 
open boundary conditions, and the QISM is still availableil. A new ^-operator is introduced 
corresponding to the boundary impurity. In most works, the i^-operator is a 2 x 2 matrix 
with the elements being c-number describing the ^-coupling to the applied fieldfi The 
Sklyanin's formalism can be extended to the generic representations of i^-operator, which 
is written as a 2 x 2 matrix but with elements being operators rather than c-numbers0@. 
This operator-valued K plays an useful role in constructing the boundary problem where 



the quantum degrees of freedom of the boundary enter interactions. Of course generally, 
both R, ^-vertices could be interpreted as the inhomogeneous vertices in a 2D lattice model. 
Our model corresponds to a very special one that the only "inhomogeneity" comes from the 
boundary row, leaving others uniform. It is built from 

K(X) = ss , R A0 (X - ic)ss>R M (\ + ic) (10) 

and 55'i?(A) A0 is given by 

S+S> I \ — k S+S ' V - V 

SS ,*-(A)=- z n j-j n (id 

l=\S-S'\ fc=|5-5'|+l n^l,n=\S-S'\ yi yn 

It is straightforward to show that the doubled monodromy matrix 

0(A) = T{X)K{X)T~ 1 (-X) (12) 

satisfies the reflection YB relations and its trace 6*(A) = tr4©(A) satisfies [9(X),9(fi)] = 

0, VA,/i. Similarly, because H oc ^ ln^(A)|A=o> the Hamiltonian (7) is indeed integrable. 

Its spectrum is uniquely determined by the following Bethe ansatz equations (BAE): 

Xj + i(S' + c) Xj + i(S' - c) Xj + iS 2N+1 = ^ A 3 - - X t + i Xj + Xi + i 

Xj-i(S' + c)' Xj-iiS' -^^Xj-iS 1 l^Xj-Xi-i ' Xj + ^-i' [ } 

The energy of the Hamiltonian (7) is 

M q 

up to a rapidity-independent constant. The magnetization is given by S z = NS + S' — M 
with M being the number of down-spins. It is recalled that when S' = S, c = (or J = J ), 
the impurity is just the one more site of the chain. For convenience, we put J = 1 in the 
following text. 

III. GROUND STATE, BOUNDARY CORRELATOR AND BOUNDARY 

STRINGS 



Due to the reflection symmetry of the model and its Bethe ansatz equation, there is a 
restriction on the rapidities: Xj ^ ±A/, for j ^ I. Therefore, Aj = is forbidden in this 
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system. Generally, the bulk solutions of (13) can be described by the following strings in 
the thermodynamic limit 

A? l7 = A? + i(n-2j + l), j = l,2,..-,n (15) 

with A" a positive real number. Since c and — c give the same Hamiltonian, not lossing 
generality, we consider only c > (c real) or Im c > (c imaginary) cases. For c < S' 
and or imaginary c, (15) are the only possible solutions of the Bethe ansatz equation (13). 
For each class of states classified by n-strings, we introduce the usual density distribution 
function p n (A) and p n ,h{X), representing occupied states (particles) and missing states (holes) 
respectively. The Bethe ansatz equation of the n-strings reads: 

OO 1 

p n , h (X) + J2 An^(A) = a n , 2S (A) + — + C* S (A)], (16) 

where a n (A) = n/27r(A 2 + n 2 /4), a n> /(A) = £™"( n,J ) a„ + i + i_ 2 A:(A); A ni is an integral operator 
with the kernel 

min(n,l) — 1 

A»,«(A) = 0| n +/|(A) + 2 ^ a„ + i_ 2 fe(A) + a| n _;|(A), 

fe=i 

0^ mp (A) = a ni 2S"(A — ic) + a n,2S"(A + ic) is the impurity contribution; <^ d9 (A) = a n (A) + 
On+i(A) +a n _i(A)(l — 5i, n ) is the surface or edge term, which is independent of the magnetic 
impurity. Notice S(X) in (p^f 9 is included to cancel the A™ = term, which is the solution 
of the Bethe ansatz equation but corresponds to a zero wave function (a direct result of the 
restriction Xj ^ A;). In the ground state, only 25'-strings exist@ and (16) is reduced to 

A 2S ,2SP2 S (X) = a 25 , 25 (A) + ^[&{X) + 2 t(A)], (17) 

By Fourier transforming (17), we readily obtain 

p 2s (X) = pl s (X) + ^[pTs P W + tifW], (18) 

*< A > - jJ^ry < 19 > 

tovrw- 1 f sin H s>UJ ) cosh(cw) _ iujX 
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^^'=^/ ta ti!h^r i) e "'^ + ^ (A) - <22 » 

The ground state energy can be readily calculated as 

Eq/N = f° ulk + jr{E? mp + F® dg + -fb ulk ), (23) 

fLi k = lm\)-n\ + s)], (24) 

FL P = \ £ + ^ - 5'| + zrc) - + 1(5 + S') + ire)], (25) 

H-^ + ^-lD-^ + ^^+l))]- (26) 

where \I> is the digamma function. As a byproduct, the correlator of 5' and 5*1 can be exactly 
derived for the present model. When 5 = 1/2 or 5' = 1/2 we have 

<S'-S X >= -^jE . (27) 

Since c (and therefore J) is only included in p^ p (\), we need only the impurity energy F® . 
The boundary correlator can be calculated as 

< S' • Si >= J~ -Q^F imp . (28) 

Now we turn to c > 5' case. In addition to the n-string solutions (15), an imaginary 
mode A = i{c — 5') appears to be a solution of the Bethe ansatz equation (13). In fact, in 
this case a so-called boundary n — fc-stringi! is possible solution of the BAE 

Ks~ k ' m = *( c ~ + m = k, k + 1, ■ ■ ■ ,n, (29) 

where k < S' — c or k = 0. Generally, there is no restriction for n in the spin chain with a 
boundary field^. However, in our case, A = ±i(S' + c) are not solutions as we can see from 
the Bethe ansatz equation. That means n < 25". In addition, for 2c = integer case, k must 
be zero due to the restriction Xj ^ ^ I. Suppose there is an n — k boundary string in 



the ground state configuration. The change of the 25'-string distribution pb s (X)/2N can be 
readily derived from the following equation 

n 

A 2 s,2sPbsW = ~ Efe^A - z(c -S' + I)) + a 25 , 25 [A + z(c - S' + I)]}. (30) 

l=k 

The energy carried by the boundary string is 

7T f n S 

ebs = -2J °WAK(A) - E g 2 _ (c _ y + /)2 - ( 31 ) 

By solving (30) via Fourier transformation and submitting pb s into (31), we find e^s = 0. That 
means the boundary string contributes nothing to the total energy in the thermodynamic 
limit. It corresponds to a charged vacuum in the sine-Gordon theory!!. We remark that 
some boundary string will be stablized with a finite magnetization. Such a situation is much 
like those of the fermion systems with boundary potentiall! or Kondo impurityBEl. 



IV. THERMODYNAMICS 

In this section, we consider the thermodynamics of c < 5" (antiferromagnetic) case. 
The thermodynamic Bethe ansatz equations can be derived by following the standard 
methodiliE. At finite temperatures, the solutions of BAE are described by (15) . The 
energy of the system takes the form 

E 

= —71 

71=1 " n=l 

where H is the external magnetic field. The entropy of the system reads 



TP OO „ OO „ 

]y = E J «n, 2S (A)p„(A)rfA + J2 nH J PnWdX, (32) 



oo „ 

S/N=Y, {{pn + Pn,h)^{pn + Pn,h) - Pn^Pn- Pn,h^Pn,h}dX- (33) 
n=l J 

By minimizing the free energy F = E — TS we readily obtain the following equation 

ln(l + Vn) = UH ~™ n ' 2S + £ A n>l ln(l + C 1 ), (34) 
1 i=i 

where 77 n (A) = Pn,h(A)/p n (A) and 770(A) = 0. With the identities 
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G, Bi >m — GB 2i m — <5i,m.G, (36) 

where B nm and G are integral operators with the kernels a n>m (A) and l/2cosh(7rA) respec- 
tively, (34) can be reduced to 

zj cosn(7rAj 



with the boundary condition 



lim ^ = * (38) 



(37) is almost the same equation to that of the 2S"-channel Kondo problem.il with only a 
different driving term. The free energy reads 

F/N = fbuik + j^F imp + —F e d g , (39) 
hulk = fU -Tj[2 cosh(TrA)]- 1 ln[l + r) 2S (X)]d\, (40) 
F edg = ±H-~Tj[2 cosh(7rA)]- 1 {ln[l + m (X)} + ln[l + 7fe(A)]}dA, (41) 

1 1 oo - 

F imp = ~(S' $r (A) ln[l + rj~ (X))d\, (42) 

Z Z n=l J 

It is not easy to reduce further F imp for arbitrary real c. The boundary behaves always as a 
spin-1/4 (in fact one half of a spin-1/2) and its critical effect in the XX Z spin-1/2 chain has 
been discussed in a previous workH. In our case, the "boundary spin" shows overscreened 
critical behavior as will be discussed in the following text. 

A. High Temperature Limit 

When T — > oo, the driving term tends to zero. Therefore in this limit the functions r) n (X) 
tend to constants rj~ satisfing the algebraic equatiorS! 

vf = + + (43) 

with the boundary conditions 
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Vq = o, lim = ^ = 2x . (44) 



The solutions of (43) are 



= sinhVj- 1V„ _ L 
smn 



Substituting (45) into (41) we get 

Fedg ~ -iTln(2cosha: ). (46) 

Notice that only the pure boundary term is given in the above expression. Obviously, the 
open boundary's contribution to the free energy is one half of a spin-1/2. The entropy of 
the open boundary is In \/2. To calculate the impurity free energy, we note that the integral 
kernel $^ mp in (42) can be written with real variable as 

n 

$r P (A) = E E «|n+i+25' + 2r C -2 fc |(A)e(n + 1 + 2S' + 2rc - 2k), (47) 

r=± k=l 

where e(x) = sign(x) and e(0) = 0. Since all a„(A) are equivalent to 5(\) in T — > oo limit, 
we can replace (47) by 

$r P (A) - a n , 25 '- C/ (A) + a n , 25 ' +C/ (A), for Cl = 2c (48) 
$r p (A) (A) + (A) + 5(A)E 0n,25'-c / +2Z-l, for a^2c, (49) 

where cj is the integer part of 2c. Therefore, the impurity free energy reads 



1 

4 

1 

4 

1 c ' 



i^mp 7 r[ln(l + 77 2S ,_ CJ ) + ln(l + ?7 25 , +C/ )], / or c ' = 2c > ( 50 ) 

*imp ~ -7 T [ lll ( 1 + V2S'-cj) + ln (! + V2S'+cj)} 



+ 9 T E[ ln ( 1 + %5'- Cj +2i-l) - ln %S'^ I+ 2I-l]. M C/^2C. (51) 
Z Z=l 

The entropy of the impurity is 

S zmp = \n^/(2S' + l) 2 -c 2 I , for Cl = 2c, (52) 
5 imp = ln^/(2S" + l) 2 -^ 

-f> ; 2S'-Cj + 2l f()r ^_^ 2c 

2=1 ^(25" - C/ + 21 - 1) (25' - a + 21 + 1) 
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^From (52-53) we see that even at high temperature limit, the impurity spin is split into 
ghost spins via the boundary bond deformation. Though the magnetization of the impurity 
is almost the same to that of a free spin S', the entropy is strongly interaction-dependent. 



B. Low Temperature Limit 

When T — > 0, the driving term in (37) diverges. That means r]2s —* and all other r] n 
tend to constant rj^ which satisfy the same equation of r]~ but with a different boundary 
condition. Since the equation is decoupled at n = 2S, we have different solutions! for n > 2S 
and n <2S 

+ sinh 2 (n - 2S + l)x 

Vn = T-T2 !> f° r n > 2S > ( 54 ) 



sinh xq 



s i n 2 E n + 1 

ti= ■ 2 2S , +1 -1, for n<2S. (55) 

S111 2(5+1) 

The residual entropy of the open boundary is 

1 7T 

Sedg = ~ hl[2 COS 2 , S+1 x \- ( 56 ) 

For 2c = C/, the impurity behaves still as two ghost spins S' + ci/2 and 5" — C//2. The 
entropy of a ghost spin S reads 

S ghost = Un[2(S-S) + l], for \S\>S, (57) 
i s i n zl M+l 

S ghost =-ln . 2 j +1 , for \S\<S. (58) 

1 S111 2(S+1) 

The summation of the two ghost spins' entropy gives that of the whole impurity. When 
cj 7^ 2c, the difference between the residual entropies for c/ ^ 2c and cj = 2c, i.e., AS imp 
reads 

1 CI 

&S imp = -- 23[ln(l + / 2 5'- Cj +2/-i) - ln/ 2 s'_ cx+2 j-i], (59) 

2 i=i 

with / n = lim a . ^ ■ This result shows that the spin configuration of the ground state is 
very complicated and strongly depends on the impurity-bulk coupling. It would be plausible 
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to coin it as a local spin glass. In fact, the residual entropy has jumps at c = cj/2. That 
means quantum phase transition occurs for c across a half integer or an integer. 

To obtain the leading order of some thermodynamic quantities such as the specific heat 
and the susceptibility, we need the low temperature (T << T&) expansion. This can be done 
by following the standard method developed for the multi-channel Kondo problem!. For 
T — > 0, only the excitations near the Fermi surface (A — > ±00) are important. The driving 
term in (37) can be approximately replaced by — (it /T)exp{— 7r|A|). We introduce the new 
variables (± = ±7rA + ln(7r/T), then rj n take the following asymptotic formi 

Vn(C±) ^Vn + ian + PnxDe-^, for n>2S, (60) 
»7n(C±)~^ + (an + ^)e~ TC± , for n<2S. (61) 

Here a n and (3 n are constants, r = 2/(5 + 1) and ± denotes the two Fermi points. For 
imaginary c = ib, the free energy of the impurity reads 

ft. ~n*-T /[ 2cosh(c+ U_ lnf) + 2cosh(c _U-m f ) ' ^ + ■»•>«■ < 62 > 

Notice that we have replaced £± by ( in the integral. In this case, the bond deformation 
does not change the effective strength of the impurity but the energy scale (Kondo 
temperature )0 

Tfc ~ 7r cosh _1 (7r6). (63) 

The system behaves as a 25-channel Kondo system with an impurity S'. For real c and 
cj = 2c, the free energy of the impurity can be rewritten as 

mp tmp 2cosh(C-lnf) 1 ' 



For cj ^ 2c, 



F imp ~F? mp -Tj2j G ± [±(( - In £)] ln(l + ms , ±CI )d( 

•^§/ [I( C -lnf)7+(c-l Cj )=» ln(1 + ^- Cl+ v-i)d(, (65) 
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where 

Notice G±(X) are convergent in the real axis since c — C//2 < 1/2. The specific heat and 
the susceptibility can be easily derived from the free energy. With different values of c, 
different quantum critical behavior may appear, (i) For S ± cj/2 > S, both the ghost spins 
are underscreened and the leading terms in the specific heat and the susceptibility are the 
Schottky term and the Curie term respectively, (ii) For S' — ci/2 < S < S" + c//2, no matter 
how large S' is, 

C imp ~T\ Ximp ~T- 1 + 0(T T - 1 ), (67) 

which indicate a novel critical behavior, (iii) For S' ± cj/2 < S, the system behaves as a 
conventional overscreened Kondo system. The above results show that the bond deformation 
has two effects: The half integer part of c (cj/2) renormalizes the impurity spin and the rest 
(c — cj/2) renormalizes the effective energy scale Tk (Kondo temperature). 



V. CONCLUDING REMARKS 

In conclusion, we propose an integrable model of a boundary impurity spin S' coupled 
with an open Takhatajian-Babujian spin-S* chain. The relation between the present model 
and the bulk impurity problem in a spin chain is discussed. In our model, when S or 5" is 
one half, The "fine tuned" effect in the periodic integrable models!'! is overcome and the 
interaction term in our case takes a very simple form. The coupling constant J can take 
arbitrary value without destroying the integrability of the Hamiltonian, while in the periodic 
models, there is a constraint to J. Though a similar s'sR A0 (^ — c) can be introduced in the 
periodic mo delJi, the parameter c must be real (imaginary in our case) which describes 
a weak-linked impurity to the bulk. The interaction only affects the energy scale (Kondo 
temperature) but does not change the fixed point of the system. With an imaginary c, the 
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model Hamiltonians constructed for bulk impurities are non-Hermitian and their spectra 
generally lie in the complex plane!! rather than in the real axis completely. In our model, 
both real c and imaginary c define a Hermitian Hamiltonian due to the reflection symmetry, 
and the coupling constant J meets all physical situations. Some new quantum critical 
phenomena driven by the impurity-bulk coupling have been found, which can never appear 
in the periodic models: (i)The stronger coupling J may split the impurity spin into effective 
"ghost spins" S' — cj/2 and S' + Cf/2. That means the coupling not only change the energy 
scales (Kondo temperature) as in the conventional Kondo problem but also renormalizes 
the effective strength of the impurity spin. It seems that the strength of these ghost spins 
does not change via temperature. Such a phenomenon reveals a pure correlation effect. We 
remark a similar effect can be induced by the impurity potential in the Luttinger-Kondo 
problems0H.(ii) Depending on the strength of the coupling, the system may show a variety of 
critical behavior differing from those of the conventional Kondo problems. A typical example 
is that when S' — c//2 < S < S' + Cj/2, the leading term in the susceptibility is Curie type, 
while that of the specific heat is overscreened 2S'-channel Kondo type. Such a fascinating 
non- Fermi liquid behavior has never been found in the conventional impurity problem (notice 
these are induced by the same impurity). (iii)The open boundary, which can be produced 
by either an impurity (no matter magnetic or non-magnetic) or bond deformation, shows 
overscreened multi-channel behavior as long as the bulk spin S > 1/2. Such an effect is 
caused by the self-avoiding of the scattering of each spin- wave with its reflection counterpart 
and is common to the multi-channel systems in ID. Our results strongly suggests that some 
new intermediate fixed point may exist for the Kondo problem in a strongly correlated host. 
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